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Abstract. - In this Letter we investigate the mechanism for overcharging of a single spherical 
colloid in the presence of aqueous salts within the framework of the primitive model by molecular 
dynamics (MD) simulations as well as integral-equation theory. We find that the occurrence 
and strength of overcharging strongly depends on the salt-ion size, and the available volume 
in the fluid. To understand the role of the excluded volume of the microions, we first consider 
an uncharged system. For a fixed bulk concentration we find that upon increasing the fluid 
particle size one strongly increases the local concentration nearby the colloidal surface and that 
the particles become laterally ordered. For a charged system the first surface layer is built 
up predominantly by strongly correlated counterions. We argue that this a key mechanism to 
produce overcharging with a low electrostatic coupling, and as a more practical consequence, 
to account for charge inversion with monovalent aqueous salt ions. 



Overcharging, or charge inversion, is defined as the situation where a charged colloid 
(macroion) accumulates close to its surface more counterions than necessary to compensate its 
own bare charge. This effect was already discovered in the beginning of 80's, both by computer 
simulations |Q] and analytical studies . Based on reversed electrophoretic mobility, some 
experimental and numerical (MD) |^ studies provide some hints for the manifestation 
of overcharging and its possible experimental relevance. More recently, it has regained a 
considerable attention on the theoretical side [pHlSl- 
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Overcharging is rather weU understood for counterions in salt- free or low sah environment, 
when excluded volume effects play no role. The underlying physics at zero temperature of 
such non-neutral systems can be quantitatively explained with Wigner crystal (WC) theory 



11 1 . The basic concept is that the counterions form a two-dimensional lattice on the 



macroion surface, and when overcharging counterions are present on this layer the energy 
of the system is lowered compared to the neutral case. This feature can be directly and 
exactly computed for a small number of counterions at zero temperature, and was illustrated 
in Refs. This WC approach remains qualitatively correct for finite temperature as long 
as the Coulomb coupling is very high. The crystal then melts into a strongly correlated liquid, 
where the local order is still strong enough to lower the free energy for the overcharged state. 

The situation becomes much more complicated for aqueous systems, where the coupling is 
weak and in addition salt is present. One-particle inhomogeneous integral equation theories 
can describe some of these situations fairly well H, |,14, 1^ but the computed correlation 
functions do not necessarily give direct insight into the physical mechanisms behind these 
effects. The presence of excluded volume interactions can lead to layering effects near the 
macroion, which are known from simple fluid theories. Here, due to the presence of charge 
carriers of both signs, this can even lead to layers of oscillating charge inversions |||,14, 1^. This 
overcharging was also observed (by integral-equation and simulation ^ ) for monovalent 
salt-ions of large size. However, until today, the basic mechanism of charge inversion for dense 
salt solutions remains unclear. 

The goal of this Letter is (i) to study in detail the role of the excluded volume contribution 
for the overcharging of a colloidal macroion in the presence of salts and (ii) to provide a 
qualitative insight into the mechanism behind these effects. We find that in the presence 
of salt the contribution of excluded volume can be so important that the size of the small 
ions dominates the occurrence of overcharging and the overcharging strength increases with 
increasing ionic size when the electrolyte concentration is fixed. Even for monovalent ions with 
high enough ionic size, we observe overcharging, which can not be explained with a salt-free 
WC picture due to the low Coulomb coupling strength. In order to have the simplest system 
we study only the cases where the colon and the counterions have the same size. This will 
reduce the effects of depletion forces which lead to nontrivial features already in neutral hard 
sphere fluids. 

Our proposed mechanism will rest on the following arguments: For a fixed salt concentra- 
tion, the available volume in the fluid is function of the electrolyte particle size. More precisely, 
the entropy of the solution is decreased by enlarging the size of the salt-ions jl^ , which en- 
hances inter-particle correlations. On the other hand the interface provided by the macroion 
leads to an increase of microion density close to the macroion, and promotes there also lateral 
ordering, even in the absence of strong electrostatic coupling, similar to a prefreezing phe- 
nomenon. However entropy alone can never lead to overcharging, since in this limiting case 
colons and counterions have the same radial distribution. But ordering and weak electrostatic 
correlations can lead to overcharging, as we are going to prove in this Letter. 

In this work we carried out MD simulations and HNC/MSA (hypernetted chain/mean 
spherical approximation) integral-equation to study the overcharging in spherical colloidal 
systems within the primitive model. In particular, the solvent enters the model only by its 
dielectric constant and its discrete structure |l9) is ignored. The system is made up of (i) a 
large macroion with a bare central charge Q — —ZmS (with Zm > 0) and (ii) symmetric salt 
ions of diameter a and valence Zc- The system is globally electrically neutral. 

For the simulation procedure, all these ions are confined in an impermeable cell of radius 
R and the macroion is held fixed at the cell center. In order to simulate a canonical ensemble, 
the Langevin thermostat has been used to predict the ions trajectory as has been done in 
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Table I - Simulation run parameters for the charged fluid (A — F) and the neutral fluid (G — H). 



parameters salt valence Zm o-/Ib f 



run A 


2:2 


48 


1 


2.3 X 10"' 


run B 


2:2 


48 


0.5 


2.9 X 10-2 


run C 


2:2 


48 


0.25 


3.6 X 10"^ 


run D 


1:1 


10 


1 


2.3 X 10"^ 


run E 


1:1 


48 


1 


2.3 X 10"^ 


run F 


1:1 


48 


0.5 


2.9 X 10~2 


run G 






1 


2.3 X 10"^ 


run H 






0.5 


2.9 X 10"2 



Refs. 0. 

Excluded volume interactions are taken into account with a purely repulsive (6-12) Lennard- 
Jones (LJ) potential characterized by the standard (cr, e) length and energy parameters, re- 
spectively, which is cut-off at the minimum [po| . The macroion-counterion distance of closest 
approach is defined as a. Energy units in our simulations are fixed by e = ksT (with T = 298 

K). _ ^ 

The pair electrostatic interaction between any pair ij, where i and j denote either a 
macroion or a microion, reads Ucoui{r) — ksTlB^^-^ , where Ib — e^/(47reoerfcs2o) is the 
Bjerrum length fixing the length unit. Being interested in aqueous solutions we choose the 
relative permittivity = 78.5, corresponding to Ib = 7.14 A. The time step is At = O.Olr 
with T = r^^, where F is the damping constant of the thermostat. Typical simulation runs 
consisted of 2 — 7 x 10^ MD steps after equilibration. 

In order to give prominence to the effect of ionic size, we choose to work at fixed distance 
of closest approach (between the centers of macroion and the salt-ion) a = 21b = 14.28 A. 
In this way, the electrostatic correlation induced by the colloid remains the same (i.e. fixed 
macroion electric field at contact) no matter what the ionic size is. Thus by changing the ionic 
size (T one also changes that colloidal radius accordingly so that a remains constant. The salt 
concentration p is given by N-/V where V = ^t^B? is the cell volume and iV_ is the number 
of colons. We will restrict the present study io p = IM salt concentration. The fluid volume 
fraction / is defined as {N++N^) (^)'^ where iV+ is the number of counterions. To avoid size 
effect induced by the simulation cell we choose a sufficiently large radius R — 8.21b yielding 
to more than 1000 mobile charges. Simulation run parameters are gathered in Table |. The 
HNC/MSA calculations were performed using the technique presented in |l^ and references 
therein. Here it is assumed that the system size is infinite, and the bulk salt concentration is 
fixed. In practice there should be no observable difference in the correlation functions (between 
HNC/MSA and MD) close to the macroions, because the wall effects die off sufficiently fast. 

We first illustrate the excluded volume correlations present for a neutral hard sphere fluid 
(runs G— H) [identical to systems A — B and E ~ F, but uncharged (see Table |)]. Although 
these kind of systems are "simple" fluids |Q , it is fruitful to elucidate what exactly happens 
at this "low" fluid density upon varying the fluid particle size a in the presence of a single 
large spherical particle. To characterize the fluid structure we consider the pair distribution 
function g{r) between the colloid and the fluid particle, which is just proportional to the local 
density n{r): n{r) = pg{r). Results are depicted in Fig. (|l|). For the large value a = Ib, one 
clearly observes a relative high local concentration as well as a short range ordering nearby 
the colloidal surface. When the particle size is reduced by a factor 2 (holding a fixed), the 
behavior is qualitatively different and the system is basically uncorrelated. By increasing a 
at fixed fluid density p the mean surface-surface distance between particles is reduced which 
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Fig. 1 - Pair distribution function g{r) between the colloid and the fluid particle in the neutral state 
(runs G — H) for two particle sizes a. The origin of the abscissa is taken at the distance of closest 
approach a = 21b- Lines and symbols correspond to HNC theory and simulation respectively. 



in turn leads to higher collision probability and thus to higher correlations. In other words, 
by reducing the available volume one promotes ordering [p^,p2[. This seems to be trivial in 
the bulk, but the presence of the colloidal surface induces an even enhanced ordering and the 
system can prefreeze (order) close to the colloid. Note that the same effects are naturally 
present for a fluid close to a planar wall pT|,p^. 

To characterize the overcharging we introduce the fluid integrated charge Q{r) which 
corresponds to the total net charge in the fluid (omitting the macroion bare charge Zm) 
within a distance r from the distance of closest approach a. Results are sketched in Fig. ||(a) 
and Fig. |^(b) for divalent and monovalent salt ions respectively. Theoretical and numerical 
analysis are in very good qualitative agreement. 

For the divalent electrolyte solutions (runs A-C) we observe that overcharging is strongly 
dictated by the ionic size a [see Fig. ||(a)]. For small ions (run C) no overcharging occurs (i. 
e. Q{r)/Zm < 1), which is a non trivial effect probably related to the forming of ionic pairs 
(counterion-coion pairing) due to the strong pair interaction of SkgT . This delicate point will 
be addressed in a future study, and was also observed in [|l5|. Upon increasing cr, the degree of 
overcharging increases. We carefully checked that the distance r = r*, where Q{r*) assumes 
its maximal value, corresponds to a zone of coion depletion (in average there are less than 2 
colons within r*). This implies that also the absolute number of counterions at the vicinity of 
the macroion surface increases with increasing a. This is qualitatively in agreement with what 
we observed above for neutral systems. However we are going to show later that electrostatic 
correlations are also concomitantly responsible of this extra counterion population (increasing 
with (t) in the vicinity of the colloidal surface. 

For the monovalent electrolyte solution [see Fig. ^(b)] overcharging occurs for a — Ib 
(runs D and E). In respect to the salt- free WC picture this is rather unexpected since here 
the "plasma" parameter Tec = IbZ^/occ-, where = (nc)^^^^ (c standing for the two- 
dimensional surface counterion concentration) is the mean distance between counterions on 
the surface, is small. More precisely for Z^ = 10 (run D) we find F^c « 0.8, and for Z^, = 48 
(run E) we have T^c < 1.0 [Q. But following this salt-free approach, it is necessary to have 
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Fig. 2 - Reduced fluid integrated charge Q{r)/Zm as a function of distance r for for three different 
particle sizes a. (a) divalent salt ions (runs A-C), (b) monovalent salt ions (runs D-F). The origin of 
the abscissa is taken at the distance of closest approach a — 21b 



at least Fee > 2 to get overcharge [g^. Note that from run D {Zm = 10) to run E {Z„i — 48) 
one increases the macroion surface-charge density leading (for Zm = 48) to a higher absolute 
overcharging Q{r*) but a weaker ratio Q{r*)/Zjn, which is qualitatively in accord with the 
WC picture, since the maximal overcharging is proportional to -y/c. A closer look on Fig. 
^(b) reveals that for Z^ = 48 (run E) r* /a is shifted to the right compared to the divalent 
case. This is merely a packing effect and it is due to the fact that for monovalent counterions 
the macroion charge (over)compensation involves twice more number of counterions particles 
than in the divalent case. Therefore for the macroion charge density under consideration 
{Zm = 48), more than one counterion-layer is needed to compensate the macroion charge. 
Again for a smaller ionic salt size a = 0.51b (run F) the overcharging effect is canceled. 

Recently a depletion of salt-ions {total local density of co and counterions) near the 
macroion surface was reported for bulk salt-ion concentrations similar to ours but very low 
surface charge density |26|. For all of our investigated cases (runs A-F), however, we never 
observe such a depletion. The reason is that such an effect vanishes as soon as the ionic size is 
sufficiently large and/or the surface-charge density is large enough, which is in agreement with 
our findings. Indeed, similarly to what happened with image charges, the self-image repulsion 
is only relevant for low charge density . Note that the parameters leading to overcharging 
cannot give rise to an observation of salt-ion depletion. 

For the WC picture to be effective we need strong lateral correlations which cannot come 
from pure electrostatic effects. To see if such correlations are present we consider the local two- 
dimensional surface counterion structure. We analyzed in our simulations the two-dimensional 
counterion pair distribution gcc{s), where s is the arc length on the macroion sphere of radius 
a. All counterions lying at a distance r < a + 0.5cr from the macroion center are radially 
projected to the contact sphere of radius a. Predominantly counterions are present in the 
first layer. For the neutral system (run G) we analyzed the small neutral species. Results 
are given in Fig. ^. We observe that all systems with an ion diameter a = Ib show their 
first peak at roughly Itr, and a weaker second peak at about 2(7, exhibiting long range surface 
correlations. The second peak is very weak in the neutral system G, but clearly visible in 
system A and E. Due to the stronger electrostatic repulsion the second peak for the divalent 
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Fig. 3 - Two-dimensional surface counterion correlation functions gcc{s) for divalent salt (run A - C), 
monovalent salt (run E), and a neutral system (run G). 



system A is somewhat further apart than for the monovalent system E. The systems B and C 
with smaller ion diameters show a correlation hole of size « 2cr, which is of purely electrostatic 
origin. Therefore lateral correlations can be produced either by pure entropy effects (run G) 
or pure electrostatic effects (run C), or we can have an enhanced lateral ordering due to the 
interplay of both (run A, E) . Qualitatively, WC arguments are still applicable to monovalent 
systems (such as run E) if one considers an effective low local surface temperature stemming 
from the strong surface ordering. 

In summary, the observed overcharging for low Coulomb coupling can be traced back to the 
complicated interplay of entropy and energy by the following two effects. First, by enlarging 
the excluded volume of the salt ions at fixed concentration, one lowers the accessible volume 
to the fluid particles. Second, the presence of the large macroion provides an interface near 
which the density of the fluid is increased compared to the bulk, and the solution can prefreeze 
due to entropic effects, provided the available volume gets low enough. The closest layer to the 
interface possesses already strong lateral correlations, even for a neutral system. If the system 
is additionally charged, then even weak Coulomb correlations can lead to the formation of 
a strongly correlated liquid, where the overcharged state is energetically favorable, as shown 
in 1^, |ll|, ^ |l^ . The order of this counterion layer is however not created by electrostatic 
interactions as in the normal WC picture, but it is largely due to entropic effects. The observed 
overcharging effect might have implications for the stability of colloidal suspensions. Additions 
of monovalent salt will eventually make colloidal suspensions unstable due to the onset of the 
van-der-Waals attractions. Upon addition of even more monovalent salt there is the possibility 
of salting the precipitate again in, as has been seen for polyelectrolyte systems The 
observation of such a reentrant transition could be an important hint towards the existence 
of overcharging with monovalent salt ions. 



* * * 
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